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Abstract 

Quantum fields written on noncommutative spacetime (Groenewold - Moyal plane) obey twisted commu- 
tation relations. In this paper we show that these twisted commutation relations result in Hanbury-Brown 
Twiss (HBT) correlations that are distinct from that for ordinary bosonic or fermionic fields, and hence can 
provide us useful information about underlying noncommutative nature of spacetime. The deviation from usual 
bosonic/fermionic statistics becomes pronounced at high energies, suggesting that a natural place is to look at 
Ultra High Energy Cosmic Rays (UHECRs). Since the HBT correlations are sensitive only to the statistics of 
the particles, observations done with UHECRs are capable of providing unambiguous signatures of noncommu- 
tativity, without any detailed knowledge of the mechanism and source of origin of UHECRs. 



Cosmic Rays with energies around 10 18 eV and higher are called as Ultra High Energy Cosmic Rays (UHE- 
CRs) 13. They are the highest energy particles known to us and often have energies 10 7 times more than that 
produced by LHC. Inspite of recent advancements (both theoretical and experimental), UHECRs pose a consid- 
erable theoretical challenge: the source and mechanism of origin of such high energy particles 4 , 5 and their 
composition are areas of active research [BHH]- 

Due to their extremely high energies, UHECRs are not only an excellent arena for testing the validity of known 
laws of physics [MTU] but are also some of the best places to look for signatures, if any, of new physics e.g. theories 
with Lorentz violation and/or deformed dispersion relations |HH15j . In this paper we aim to show that UHECRs 
can be used to look for signatures of a particular model of nonlocalities coming from the underlying noncommuta- 
tive nature of spacetime at short distances. 

Simple intuitive arguments involving standard quantum mechanics uncertainty relations suggest that at length 
scales close to Plank length, strong gravity effects will limit the spatial as well as temporal resolution beyond some 
fundamental length scale (L w Planck Length), leading to space - space as well as space - time uncertainties [16] . 
One cannot probe spacetime with a resolution below this scale i.e. spacetime becomes fuzzy below this scale, result- 
ing into noncommutative spacetime. Hence it becomes important and interesting to study in detail the structure 
of such a noncommutative spacetime and the properties of quantum fields written on it, because it not only helps 
us improve our understanding of the Planck scale physics but also helps in bridging standard particle physics with 
physics at Planck scale . 

This noncommutative scale need not always be as small as Planck scale. For instance if there are large extra 
dimensions, then the "effective Planck scale" can be at much lower energies (usually taken between TeV scale and 
GUT scale, depending on specific models) , resulting in appearance of noncommutativity at scales much larger than 
Planck scale. Such large scale noncommutativity is of particular interest as its effect can be detected in present 
day or near future experiments. In this paper we are interested in looking for signatures of such large scale non- 
commutativity in UHECRs. 

Of the various approaches to model the noncommutative structure of spacetime, the simplest is one where the 
coordinates satisfy commutation relations of the form 
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; fi, v = 0,1,2,3 and 9 is a real, constant, antisymmetric matrix. (1) 

The elements of the 9 matrix have the dimension of (length) 2 and set the scale for the area of the smallest 
possible localization in the fi — v plane, giving a measure for the strength of noncommutativity [17] • The algebra 
generated by x^ is usually referred to as Groenewold - Moyal (G.M) plane [IS], and in this paper we will restrict 
our attention to this noncommutative spacetime. Equivalently this noncommutative nature of spacetime can be 
taken into account by defining a new type of multiplication rule (* product) to multiply functions evaluated at 
same point i.e. 

f(x)*g(x) = f(x)ei t ^-g(x). (2) 

One particularly important feature of G.M. plane, which makes it quite suitable for writing quantum field the- 
ories on it, is the restoration of Poincare- Hopf symmetry as Hopf algebraic symmetry, by defining a new coproduct 
(twisted coproduct) for action of Poincare group elements on state vectors [TO] [ID] [21] ■ 

Twisting of the coproduct has immediate implications for the symmetries of multi-particle wave functions 
describing identical particles [22] . For example, on G.M plane the correct physical two-particle wave functions are 
twisted (anti-)symmetrized ones and are given by 

where <f) and ip are single particle wavefunctions 
operator associated with exchange given by 



= (^^)(0®VO (3) 
(of two identical particles) and Tg is the twisted statistics (flip) 



r e = T-WoT, T = ei^" 9 ^ 9 ". (4) 

Here tq is the commutative flip operator : To (<f> <8> ip) = ip ® <j). 

The above analysis can also be extended to field theories on G.M plane resulting in changed commutation 
relations between creation/annihilation operators [22], which now become 

a p o q = 77e vA9 a q a p , a p a q = r)e lpAq a\al 1 where p t * = (p°,p), p A q = p^q^, , . 

a p a q = 7?e -lpA9 a q a p + (27r) 3 2po<5 3 (p — q) J and rj = ±1 for bosons/fermions. ^ ' 

In literature there exist another approach to quantization of noncommutative fields 23 . In this approach, the 
quantization is done according to the usual rules and the quantum fields follow usual bosonic/fermionic statistics. 
However, such a quantization scheme docs not preserve the classical twisted Poincare invariance and suffers from 
UV/IR mixing 24 . In this paper we only discuss the twisted quantization as it preserves twisted Poincare invari- 
ance in noncommutative quantum field theories. 

Because of (0 the quantum fields written on G.M plane, unlike ordinary quantum fields, follow a unusual 
statistics which we call as twisted statistics. Twisted statistics are a unique feature of fields on G.M plane and 
can be used to search for signals of noncommutativity: because of the twisted commutation relations, interesting 
new effects like Pauli forbidden transitions 25, 26 can arise. The effect of twisted statistics also manifests itself in 
certain thermodynamic quantities [27j . [28) . The two-point distribution functions remain unchanged 

<a pi a P2 ) = 2( Pl ) Nl T U 3 ( Pl -p 2 ) (6) 

where = — is the thermal distribution, but for example the quantity (a pl a p2 a P3 a P4 ) gets changed 

[27]. 
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(a pl a p2ap3ap4 ) = 2 ( Pl ) 2 (p 2 ) Ng>Ng> [S 3 ( Pl - p 4 )<5 3 (p 2 - Pa) + r, S 3 (p a - p 3 )5 3 (p 2 - p 4 )] (7) 
The above differs from the commutative expression by the appearance of the factor e lpiAp2 in the second term. 
One can easily check, following a analysis similar to that done in [27], that ((6|) and ([7]) are true not only for 

(T) 

thermal distribution N p but for any arbitrary wavepacket /(p) i.e. 

< flpl a P2 ) = 2(p 1 ) /(p 1 )J 3 (p 1 -p 2 ) (8) 

and 

(a piap2 a P3 a P4 ) = 2 (px),, 2 (p 2 ) /( P i)/(p 2 ) [<5 3 ( P i - p 4 )5 3 (p 2 - Pa) + ve lpiAp2 S 3 ( Pl - p 3 )S 3 (p 2 - p 4 )] . (9) 
Here we discuss the consequences of to the HBT correlation functions. 

Hanbury-Brown Twiss (HBT) effect [29] is the interference effect between intensities measured by two detectors 
when a beam of identical particles is projected on them, with the intensities recorded by the two detectors operating 
simultaneously. This intensity increases for bosons (and decreases for fermions) when compared with the intensities 
recorded by the same two detectors, if only one is operated at a time. The correlation function for HBT effect is 
defined as 

Ji, I2 = intensities recorded by the two detectors respectively, when both are operated simultaneously. 
Zj = intensity recorded by first detector when the second detector is not operating. 
I 2 = intensity recorded by second detector when the first detector is not operating. 

The HBT correlation function C obeys 



C = 1 for distinguishable particles 
C > 1 for bosons (bunching effect) 

C < 1 for fermions (anti-bunching effect) (11) 

In the commutative case, for a beam of identical (massless, scalar 0) bosons, the HBT correlation function can 
be written as [3TJ1I3T] 

(B) (^ ) (yi)^" ) (y 2 )4 +) (y 2 )^ +) (yi)) 
' = — ^ ' (12) 

(^ _) (yi)^ +) (yi)) (^ ) (y 2 )^ +) (y 2 )) 

where yi and y 2 are the position of the two detectors, </>o + ^(y) = / ^J)3 sir c v e * P y ^ s ^ e positive frequency 

part and 4>n~\y) = J jk~ c t e%vy i s the negative frequency part of the bosonic quantum field 0o 0. Notice 
that we have deliberately taken massless fields because we are interested in looking for HBT correlation functions 
of ultra-relativistic particles. 

To account for the uncertainties in the energy measurements, which at present are quite significant for UHE- 
CRs, the incoming beam should be taken as a wavepacket, instead of plane waves. The choice of an appropriate 

1 Scalar bosons are taken for sake of simplicity but the analysis presented here can be easily generalized to higher spin bosons. 

2 We denote the usual bosonic (fermionic) fields by <j>o (tf>o) and their twisted counterparts are denoted by tf>g (ipe)- Also, the usual 
creation/annihilation operators are denoted by Cp and c p whereas the twisted ones are denoted by aj, and a p respectively. 
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wavepacket is potentially the only place where the information about the details of production mechanism or source 
of origin of UHECRs can reside. Central Limit Theorem tells us that the mean of a sufficiently large number of in- 
dependent random variables, each with finite mean and variance, will be approximately distributed like a Gaussian, 
and hence it is a good first approximation to take the wavepacket as a Gaussian wavepacket /(p) = Ne~ a ^ p ~ Po ^ 
centered around some mean momentum po- 

Taking the wavepacket to be this Gaussian, restricting ourself to only coincidence measurements and using the 
standard integrals [32] , the correlation function turns out to be 



C^ = l + e-^ (13) 

As clear from (fl~3|) . Cq B ^ depends only on y = |yi — y 2 the separation between detectors and on a. There is 
no dependence on the mean momentum po |33U34| . 

Similarly, for a beam of identical (chiral) fermions, the HBT correlation function can be written as 



(F) (^" ) (yi)7oV^ ) (y2)7oi(i±75)4 +) (y 2 )i(i±75)4 +) (yi)) 

(J\ b ) — __\ l__ HQ 

(^ ) (yi)7oi(i±75)^ +) (yi))(^ ) (y2)7oK 1 ±75)V^ +) (y 2 )) 

where, as before, yi and y 2 are the position of the two detectors, ipo + \y) = / JE~ S s u s,p c s , P e tp y is 

the positive frequency part and tp (y) = J ^J^ 3 jg- J2 S "s,p c t p e_lp y i s the negative frequency part of the 
fermionic quantum field tpo 

Since the whole analysis is done keeping ultra-relativistic particles in back of our mind, here we look for cor- 
relation between only chiral fermions (left- left or right- right): at such high energies, the particles are effectively 
massless and hence chiral fermions are the more appropriate ones to deal with H. 

Taking the incoming beam as a Gaussian wavepacket /(p) = iVe~ Q: ( p ~ Po -' centered around some mean mo- 
mentum po and restricting ourself to only coincidence measurements, the correlation function turns out to be 



a 



(F) 



9na 



- 2q p° 2 (y 2 + 4a 2 p 2 ) 1 F 1 



5 -1 



'A 2? 4^ (y Z - 4 " 2 Po 2 - 4ioy • p ) 



5 —1 

2; -5 ^ (y 2 - 4a 2 p 2 + 4^ay • p ) 



(15) 



where y = y 2 — yi is the separation between the two detectors and \F\ (a; 7; z) is the degenerate hypergeometric 
function. 



We observe that in this case, the correlation function depends on the separation between detectors y, on the 
mean momentum po of the wavepacket and the angles between y and po ■ 

In the noncommutative case, there are two important differences. Firstly, the (•) product between fields eval- 
uated at the same point has to be replaced by (*) product and secondly the expectation value is changed in 
accordance with ([7]), as the quantum fields are now composed of twisted creation/annihilation operators. Hence in 
noncommutative case for twisted (massless, scalar) bosonic particles, we have 



r ( B) = (^~>( yi ) {^-)( ya ) * y2 ^ (+) (y 2 )} * yi (+) (yQ) , 
(0 ( - ) (yi)* yi ^ +) (yi))(0 ( - ) (y2)*y 2 ^ + Hy2)) 1 ' 



3 Cq F ' = 1 between particles with opposite helicities (i.e. between left-right or right- left), as at ultra-relativistic energies, they 
behave like distinguishable particles 
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where * y = e 



In rest of the paper we restrict to only space-space noncommutativity i.e. we take = 0^ = 0. As in 
commutative case, we take the incoming beam as a Gaussian wavepacket /(p) = iVe _Q ( p_Po * ) centered around 
some mean momentum po and restrict ourselves to only coincidence measurements. C\ then turns out to be 



ci B) = 1 



4a 



4a 2 + A 2 



exp 



4a 2 y 2 - (y ■ A) 2 



2a(4a 2 + A 2 ) 



exp 



2a{p 2 A 2 - (po • A) 2 } 



4a 2 + A 2 



exp 



4ay • (p x A) 



4a 2 + A 2 



(17) 



where we have defined #y = EjjfcAfe. In getting (|17[) we have used the standard result that, for n-dim column 
matrices X and B and a n x n positive definite, symmetric square matrix A , we have 



d n X,e 



-X 1 AX+B 1 X 



( ^ 



Vdet A 



e kB T A-^B 



(18) 



We see that (|17|) not only depends on y (the separation between detectors) but also on the mean momentum 
Po, the noncommutative length scale A, as well as on the angles between A and y and between A and po- Moreover, 
as a check we can see that, in the limit A — > we get back (|13p. 

Similarly, for twisted (chiral) fermionic particles the noncommutative HBT correlation function is given by 

^ ( ' ) (yi)7o{^ (_) (y2)7o* y2 i(i±75)V> (+) (y2)}* yi i(i± 75 )^ (+) (yi) 



C. 



(F) _ 



V> ( ) (yi)7o* yi i(i±7 5 )V' (+) (yi))(^ ( ) (y 2 ho* y2 i(i± 75 )V' (+) (y2)) 



(19) 



Again restricting ourselves to only space-space noncommutativity, using Gaussian wavepackets and considering 
only coincidence measurements, C\ can be written as 



a 



(F) 



-2ap 



(K) 



2 

9na 



D -2ap 



5 zi 2 



^l 2 '2-4^ 



(Zl)a 



(20) 



where zi = y — 2iapo and Z2 = y + 2iapo- 

Expanding this in 9 and taking terms upto second order we get 



(F) 



1 



Po-(yxA) A 2 , y 2 A 2 (y-A) 2 Po 2 A 2 , (p • A) 2 , [p • (y x A)] 



121 



4a 2 



2a 



- ( —\ e- 2ap ° 2 
\9ira J 

(y 2 + 4a 2 p 2 ^ 



1 



2' 4a 



1F1 2; -; - — (y z - 4a^p ^ - 4iay • p ) iF x 2; -; -— (y^ - 4a z p z + May ■ p 



2a 
1 



2a 2 



2' 4a 



25a' 



{(4ap • (y x A) - 2A 2 ) (y 2 + 4a 2 Po 2 ) + (y ■ A) 2 + 4a 2 (p ■ A) 2 } 



7 1 



2' 4a 

(. 

175a 3 

,2 , / i„2„ L> 



1 F 1 3; -; - — (y z - 4a z p z - 4iay • p ) i-Fi 3; -; - — (y z - 4aW + May ■ p 
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2' 4a 



+ & 2X2 - 4 « 2 Po 2 A 2 - 4ia(p • y)A 2 - (y • A) 2 + 4a 2 (p ■ A) 2 + 4ia(p • A)(y • A)} 



(y 2 + 4a 2 p 2 ) 1J Fi ( 4; -; -— (y 2 - 4a 2 p 2 - May • p ) ) i-Fi ( 3; -; -— (y 2 - 4a 2 p 2 + May • p„) 



7 1 



2' 4a 
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175cr 



{y 2 A 2 - 4a 2 p 2 A 2 + 4za(p • y)A 2 - (y • A) + 4a 2 (p • A) - 4 J a( Po • A)(y • A)} 



(y 2 + 4a 2 p 2 )i J Fi f 3; ? -; ~ (y 2 - 4a 2 p 2 - May ■ p )\ 



4 We have assumed only space-space noncommutativity for calculational simplicity but one can do similar analysis with both space- 
time as well as space-space noncommutativity 
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1*1 ( 


I- 9 - 

' 2' 


1 

4a 


iFi ( 


A- 9 - 
'2' 


1 

4a 



' T ^{po-(yxA)} 2 (y 2 +4aW) 
9 1 

(y 2 - 4a 2 p 2 + 4zay • ; 



+ 0(A 3 ) (21) 



where again in writing (1211) we have defined #y = ey/sAfe- 



As can be clearly seen from (|2Tj) . the noncommutative HBT correlation function not only depends on y and po 
but also on A and various angles that A makes with y and po- 

A comparison of ([TS|) - (|T7]) and (fl~5|) - (f5T|) tells us that the HBT correlation function for twisted bosons/fcrmions 
are different from those for ordinary bosons/fcrmions and become more pronounced with increasing momenta of 
incident particles. Since A is expected to be a very small quantity, the deviations will be highly suppressed and one 
has to look at particles with high enough energy- momentum (w.r.t noncommutative scale), so that the deviations 
get sufficiently enhanced to be detectable. Hence, as stated in beginning of the paper, the best place to look for 
the signatures of twisted statistics, is to look at particles in L.H.C or in Ultra High Energy Cosmic Rays (UHE- 
CRs). Moreover, since HBT correlations are essentially quantum correlations between free identical particles and 
are manifestations of the particular statistics followed by the identical particles of a given beam, they are sensitive 
only to the statistics obeyed by these particles. Therefore despite our present lack of detailed knowledge about 
production mechanism and source of origin of UHECRs, a study of HBT correlations may still be able to provide 
unambiguous signatures, of underlying noncommutative nature of spacetime. 

The graphs shown below highlight the difference between Cq - Cg B ^ and Cq - Cg Q 




Figure 1: In the above figures, the blue and red line represent the commutative and noncommutative HBT correla- 
tion functions respectively. The figures (a), (c) are plotted for ordinary and twisted bosons and (b), (d) are plotted 
for ordinary and twisted fermions. The figures (a), (b) are plotted with |po| = 6 x 10 19 eV (i.e. at G.Z.K cutoff), 
a = 2.04 x 10 19 cV -2 (i.e taken numerically same as the present error in estimation of energy of UHECRs) [2], 
|A| = 1.47 x 10~ 24 m 2 and taken along z-axis and the angles as 9\- y = \, 8\-p a — § > 4>x-y = § an d <fi\-p = f • 
The figures (c), (d) are plotted with |y| = 1000m and same values for rest of the parameters. 

5 The persistence of correlations to large distances is attributed to the large uncertainties in our present determination of energy- 
momentum of UHECRs. With better and more precise knowledge of the energy-momentum, the correlation will decrease significantly 
resulting in better bounds on noncommutative deviations 
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Of particular interest are the various angular dependences of noncommutative correlation functions (|17p and 
(|2Tj) which are completely absent in commutative correlation functions (|13[) and (|T5|) . Due to rotation and rev- 
olution of earth, the noncommutative correlation functions will show periodic oscillations completely absent in 
commutative correlations and hence will perhaps provide the best and most unambiguous signal for underlying 
noncommutative structure of spacetime. Therefore, as claimed earlier, the information about noncommutative 
structure of spacetime, can be extracted out, from observing the nature of variation of HBT correlation functions 
(in particular by looking at angular variations) with varying certain experimentally measurable quantities, in a 
very unambiguous way. 

Also it is worth noting that if we take noncommutative length scale same as Planck scale i.e. 10 _35 m then 
due to an upper limit on momenta of UHECRs (GZK cutoff), the deviations will turn out to be O(10~ 50 ) which 
are too small to give any unambiguous signatures of it. So the noncommutative deviations are detectable only if 
the effective noncommutative scale is much larger than Planck scale, which is likely to happen in presence of large 
extra dimensions. Thus the noncommutative deviations in HBT correlations effectively provide us signatures of 
large extra dimensions. 
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